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PACS. 61.43.Hv- Fractals; macroscopic aggregates (including diffusion-limited aggregates). 
PACS. 47.53. +n- Fractals. 

PACS. 05.40.+j - Fluctuation phenomena, random processes, and Brownian motion. 



Abstract. - I analyze the topological structures generated by diffusion-limited aggregation 
(DLA), using the recently developed "branched growth model". The computed bifurcation 
number B for DLA in two dimensions is B » 4.9, in good agreement with the numerically 
obtained result of B ~ 5.2. In high dimensions, B — » 3.12; the bifurcation ratio is thus a 
decreasing function of dimensionality. This analysis also determines the scaling properties of 
the ramification matrix, which describes the hierarchy of branches. 



Nature creates an astonishing variety of branched structures. Some of these structures, 
such as trees, are created by biological systems; others, such as river networks, are created 
principally by physical phenomena. The simplest and best-understood physical model for 
the formation of branched structures is diffusion-limited aggregation (DLA) , introduced over 
fifteen years ago [fjj . The complex fractal structures generated by this model are seen in natural 
systems whose growth is controlled by diffusive processes [0 . 

Physicists have concentrated on understanding and predicting the fractal and multifractal 
properties of DLA ^ ffl. However, several authors have advanced an alternative approach, 
focussing on the topological self-similarity of DLA clusters, as measured by quantities typically 
used to describe river networks and other branched structures This Letter is devoted 

to the computation of these quantities, using the "branched growth model" which this author 
and his collaborators have exploited to compute a number of properties of DLA . Since this 
model is based on an analysis of the competition of branches in a hierarchical structure, it is 
peculiarly suited to the computation of topological quantities related to that structure. 

The growth rule for DLA can be defined inductively: introduce a random walker at a 
large distance from an n particle cluster, which sticks irreversibly at its point of first contact 
with the cluster, thereby forming the n + 1 particle cluster. Clusters grown in this way have 
an intricate branched structure, in which prominent branches screen internal regions of the 
cluster, preventing them from growing further. In addition, they contain no loops, since the 
particle closing a loop would have to attach to two pre-existing particles. The scaling of the 
radius of the cluster r with the number of particles n determines the fractal dimension D of 
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Fig. 1. - A rooted tree, showing the Strahler indices and parent Strahler indices of each branch. The 
first number is the Strahler index of a branch, and the second is the Strahler index of the branch it 
meets when it loses its identity. The total number of index 1 branches is Ni — 10, of index 2 branches 
is N% = 3, and of index 3 branches is A3 = 1. 



the cluster, n oc r . 

A number of quantities have been developed to describe the topological properties of river 
networks, which like DLA clusters have no loops. The most fundamental such quantity was 
introduced by Strahler, as a refinement of a scheme proposed by Horton Consider a rooted 
tree structure. The leaves of the structure are assigned the Strahler index i = 1. The index 
i of any other branch, which bifurcates into two branches with Strahler indices i\ and i?, is 
determined by 



Thus a change in the index of a branch occurs only when it joins another branch of roughly 



Once we have defined the Strahler index of each branch, we can also obtain N, the number 
of branches of index i. Note that a branch with several side-branches of lower index counts 
as a single branch for this numeration (fig. |l|). The bifurcation ratios Bi are then defined by 
Bi = Ni/Ni + i; by definition Bi > 2. River networks are often topologically self-similar, with 
Bi ss B and 3 < B < 5, even though they do not display conventional geometrical fractality, 
since they extend close to every point in a basin. DLA clusters in two dimensions also exhibit 
topological self-similarity, with Bi w 5.2 for i > 1 f| Q, §]. 

If the characteristic size (defined in whatever way is convenient) of branches of index i is 
Li, then we can define a length ratio Aj by A, = Lj/Lj+i. Self-similarity would then demand 
that Ai ss A, with the fractal dimension given by D = InB/ In A ||. 

The Tokunaga-Viennot ordering refines this Strahler ordering, by defining two indices i,j 
for every branch [|l0| O. The first index i is the Strahler index of the branch, while the 
second index j is the Strahler index of its "parent" (see fig. |l|). By definition, j > i. The 
ramification matrix Rij is the number of branches with the index pair ij; clearly . Rij = Ni. 
By contrast, J2i Rij gives the total number of side-branches attached to index j branches. 




h + 1, if h = h 

max(ii, if i\ ^ ii- 



(1) 



the same size (see fig. 0). 
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Since sidebranches are separated by microscopic distances on any particular branch, 

Y^Rij/NjcxLj. (2) 

i 

In addition to river networks, these quantities have also been used in computer graphics, 
and in the mathematical analysis of random binary trees jn], |l2) . These latter have B = 4 
asymptotically -we shall see below that this result is bracketed by the DLA results for d = 2 
and d —> oo. 

The branched growth model. The branched growth model for DLA is based upon a 
quantitative analysis of the competition between neighboring branches (indexed by 1, 2) that 
join at a particular node or branch point. Let the number of particles of branches 1 and 2 be 
ni, ri2, and their total growth probability be p\, p2- Let us define normalized parameters x, y 
by x — pi/ {p\ + P2) and y = ri\j{n\ +712). The total number of particles in the two branches 
is rib = n\ + ri2- We can write a simple equation of motion for y: 

dy dUl ni =x-y. (3) 



d In rib dnf, rib 

Note that the right-hand side is a function of a; and y, but not explicitly of rib- I have computed 
a similar equation of motion for x, 

<L " G(x,y), (4) 



d In rib 



by a renormalization method for two-dimensional DLA: G(x, y) is a function of x and y alone, 
but is too complicated to give here j|. 

Solving these differential equations for x, y, one generates a competition dynamics for these 
two branches in the space of their relative growth probabilities and masses, (x,y). This 
competition dynamics is symmetrical about a hyperbolic unstable fixed point at (x,y) = (^,5); 
as the number of particles rib in the two branches increases, (x, y) flows to a stable fixed point 
at either (0, 0) or (1, 1). The eigenvalue of the unstable direction at the central, unstable fixed 
point is termed v, so that for trajectories starting near this fixed point 

y(n b ) - - oc <n%, (5) 

and similarly for x. The hyperbolic nature of the unstable fixed point implies that there is a 
stable manifold flowing into it, and an unstable manifold flowing out of it and into the stable 
fixed points. 

In order to understand the statistical distribution in x — y space of branch pairs, we must 
understand not only their evolution but also their birth. I assume that branch pairs are created 
by microscopic tip-splitting events, which are indifferent to the large scale structure created 
by eqs. (pLQ). Thus those branch pairs born in the neighborhood of the unstable fixed point 
will have a roughly constant initial distribution. These branch pairs will remain "active", i.e., 
far from the stable fixed points, until ^ 1, and thus they will determine the large-scale 
structure of the cluster. Let us suppose that the distance from the unstable fixed point (or 
more properly, the stable manifold) at birth, where ^ 1, is given by e", thus defining e. A 
probability distribution for e of p(e)de = pot u ~ 1 de for e« 1 (with p a constant) implies the 
desired constant distribution near the stable manifold in the x — y plane. A simple argument 
then gives v = D^ 1 0. 

Introducing a new parameter 77 = crib, the branch competition determines a parameteriza- 
tion x(i]), y(rf) of the unstable manifold; these functions obey eqs. (^,^) with x(0) — | and 
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y(0) = |. For simplicity we choose x, y so that lim^oo x(rf), y{rj) = 0; thus for large (x, y) 
refers to the weak member of the branch pair, and (1 — x, 1 — y) to the strong member. 

We know the function G(x, y) in two dimensions, but not in higher dimensions. In this 
case, a useful approximation is given by "model Z" , for which the unstable manifold of the 
equations of motion is given by a straight line which commences at the unstable fixed point 
and continues until it reaches x — (or x = 1) |Q. It then moves parallel to the y-axis into 
the stable fixed point at either (x,y) = (0,0) or (x,y) = (1,1). Thus model Z is given by 



x(rj) = | — rf for rf < rf c = i , and eq. (|3|) implies 

,n " \v/v V>Vc w 

where fj = {l/2) l+v ~ vj{\ + v). Model Z specifies a one-parameter family of unstable 
manifolds; it is possible to determine which of these best approximates DLA as a function 
of spatial dimensionality by a self-consistency argument . 

Bifurcation numbers. - Consider a branch with a constant bifurcation number of B whose 
root has a Strahlcr index of I. We term the Strahler number of this branch to be /, the Strahlcr 
index of its root. For a DLA cluster, the total number of particles in the cluster n will be 
proportional to the total number of elementary sub-branches n e , with n e w B I ~ 1 . Inverting 
this relationship, we find 

I(n) — 1 ^ In nj In B = a In n, (7) 

where the relation is approximate because the Strahler number is by definition an integer. 
Note that we are assuming that the Strahler number of a branch is determined uniquely by 
the number of particles in the branch. In practice, the Strahler number of branches of given 
size will fluctuate in the ensemble of DLA clusters; little is known about such fluctuations. 

The main branch of the cluster is defined by starting at the root and at every branching 
choosing the stronger of the two sub- branches. This main branch eventually ends, as do all 
branches, in an elementary sub-branch tip. Equation (W) implies that 



odnn = ^pfc, (8) 



k 

where the sum is over all nodes of the main branch, and is the probability that the Strahler 
index of the main branch increases by 1 at the fc'th node, as will be the case if the two 
sub-branches at this node are close enough in size to have the same Strahler number. 
Suppose that the Strahler number of a branch is taken to be exactly 

I(n) =Int(alnn), (9) 

where Int indicates the integer part. Now consider two sub-branches of the main branch of 
size n\ or ri2- If lnri2 — lnni| > a , then they have different Strahler numbers. Otherwise, 
the probability p(n\, n?) that they have the same Strahler number is 



p(nx, na) = 1 — a\ lnri2 — lnm|. (10) 



Thus we can write 



Pfe = /y de k p(e fe )^l - a{ ln[l - y(e k n k )} - \ny(e k n k )}^j \ , (11) 
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where 3> 1 is the total number of particles below (away from the root) the fc'th node, and 
e is defined by ln[l — y(irik)] — lny(en.fc) = a -1 . The brackets (•} indicate averaging over the 
random variables at all nodes of the main branch excluding node fc, whose average is indicated 
explicitly. For nj, ;§> 1, we can use the small e result p(e) = po^ 1 to obtain 

Pk = (fy f Q dn rf- X (\ - a{ ln[l - »(„)] - In = ^ T(a), (12) 
with ln[l — y(fj)] — lny(fj) = oT x . In rcf. |l5[ ] we obtained the useful identity 

(^2^j = -{l A? Ml -»(»/)]} lnn = Alnn. (13) 

Thus, combining eqs. (|8], |l2] , |l3| ) , we obtain 

a = XT (a), (14) 

allowing solution for a. 

For DLA in two dimensions, for which yiv/) is known, solving eq. ( |l4| ) gives a fss 0.628, 
or B « 4.92, in reasonable agreement with the numerical result of B « 5.2 @, 0]. In higher 
dimensionalities, we must rely upon model Z. 

It is interesting to consider model Z in the limit v — * 0, or D = v^ 1 —* oo. In this limit, 
/ dr] "ff~ x ln[l — y(r\)\ is dominated by r\ <C r? C) and 



2i/ 

1 - ln2 + 0(i/ 2 ) 
Computing T(a) with y(rf) given by eq. @) yields 



r|a S J-)=-^J-ln 
In £7 2r/ InB 



1 - 



B- l x J 



(16) 



and eq. (|l4j) gives 



1 + Jl - exp(ln2 - 1) , , 

lira B = v yy - w 3.12 (17) 



i/-»0 



1 - ^/l - exp(ln2 - 1) 



Our somewhat surprising result is that B is a decreasing function of dimensionality. A simple 
extension of the above calculation predicts B — > oo as dimensionality d — > 1 . 

Ramification matrices. To compute ramification matrices, we must determine the 
probability that a branch of index i\ will have a side-branch of index i^. At a particular 
node, say the fc'th node, we can say approximately that 

21 — 1 21+1 

i 2 - ii = / if — < - lny fc < — . (18) 

la la 



for / > 1. Let us define rji by 



A branch with n particles will have a total length ~ n}/ D = n v . Sidebranches are separated 
by microscopic distances, thus the total number of sidebranches of a branch of size n is also 
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oc n v . The total number of branches of order i is oc B 1 , leading to 



For 77 sufficiently large, x(rj) ~ 0. To further simplify matters, we suppose that x = exactly 
for 7] > r)i, as would be the case in model Z for t]j > r\ c . Then eq. (Q) gives y{rj) — fj/f), and 

(2I-1\ . . 

7j/ = r/exp — , (21) 



2a 

which we expect to be approximately true for non-model Z trajectories. Thus 

Ri j;j -i=i oc B" 4 exp (vl/a) = B vI ~\ (22) 

This result is valid in arbitrary spatial dimensionality. In two dimensions, we obtain B v w 2.6, 
in agreement with ref. ||, which claims Rij.j^i—j oc (2.7) 1 ; however, since these latter results 
are based on a variant of the Strahler ordering, it is not clear that they are relevant. Comparing 
with eq. (Q), we see that A = B v and ln_B/lnyl = = D, as predicted ||. 
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